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We consider glass states of several disordered systems: vortices in impure superconduc- 
tors, amorphous magnets, and nematic liquid crystals in random porous media. All these 
systems can be described by the random-field or random-anisotropy O(N) model. Even 
arbitrarily weak disorder destroys long range order in the O(N) model. We demonstrate 
that at weak disorder and low temperatures quasi-long range order emerges. In quasi- 
long-range-ordered phases the correlation length is infinite and correlation functions obey 
power dependencies on the distance. In pure systems quasi-long range order is possible 
only in the lower critical dimension and only in the case of Abelian symmetry. In the 
presence of disorder this type of ordering turns out to be more common. It exists in a 
range of dimensions and is not prohibited by non- Abelian symmetries. 

1. Introduction 

Solids resist small perturbations and possess some (e.g. crystal) order. Liquids and 
gases are not ordered and offer no resistance to a static shear stress. An intermediate 
class of substances which possess ordering but respond strongly to small external 
disturbances is known as soft matter. The fourth and last possibility, no ordering 
and weak response to weak perturbations, is represented by glass phases of strongly 
disordered systems. Where should weakly disordered systems be put in this classi- 
fication? The answer depends on their symmetry with respect to transformations 
of the order parameter. If it is discrete, weak disorder is irrelevant and the system 
belongs to the class of its pure analog (e.g. the ground state of the pure and weakly 
disordered Ising ferromagnets is the same 'spin solid'). On the other hand, if the 
symmetry group is continuous then even arbitrarily weak disorder may lead to the 
formation of a glass state. However, recent studies suggest that in many cases glass 
phases of weakly disordered systems are qualitatively different from glass states at 
strong disorder and can be considered as a new type of soft condensed matter (i.e. 
combine ordering and strong response to weak disturbances). A common feature of 
these novel glass phases is quasi-long range order (QLRO). 
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As disorder is strong, the order parameter depends only on the local potential 
of impurities and hence ordering is absent. In the case of weak disorder there is a 
competition between internal interactions that tend to establish long range order 
at low temperatures and randomness that works in the opposite direction. Naively 
one could expect that stronger internal interactions win and weak disorder has no 
pronounced effect on the system. This point of view is supported by the mean 
field theory and for many years was generally accepted in the field of amorphous 
magnets. However, it is not true and even arbitrarily weak disorder may be sufficient 
to destroy long ranger order. This fact was first understood in the context of 
vortex lattices in impure superconductors^] but is valid in any system of continuous 
symmetry! The explanation of the effect is based on the low energy cost of large- 
scale collective excitations (Goldstone modes). In three dimensions this energy cost 
turns out to be lower than the energy gain due to interaction of the large-scale 
excitations with disorder. Here the continuous symmetry which is a prerequestive 
for the existence of Goldstone modes is crucial. In the systems of discrete symmetry 
weak disorder is indeed not important. Note that if disorder is weak the correlation 
length is large. This sometimes led to a wrong interpretation of experimental results 
about weakly disordered systems as evidences of long range order. 

Two questions about impure systems of continuous symmetry immediately arise: 
Whether the phase transition between high- and low-temperature phases survives 
in the presence of disorder, and if yes, what is the nature of the low-temperature 
state? Recently these questions have attracted a renewed experimental and the- 
oretical interest. On the one hand, it was stimulated by new types of disordered 
systems: superfluid Helium and liquid crystals in low-density aerogels. On the other 
hand, an important progress was achieved in studies of vortex lattices in impure 
superconductors after high T c superconductivity was discovered. In particular, it 
was found that there are two different glass states in disordered superconductors 
(Vortex glass and Bragg glass) Hta The Bragg glass state is observed at weak dis- 
order and the properties of this phase are much closer to the pure Abrikosov state 
than those of Vortex glass. It turns out that in contrast to Vortex glass, Bragg 
glass possesses topological order. Topologically ordered glass phases are predicted 
in some other weakly disordered systems and are the subject of the present review. 

Since there is no long-range order in impure systems of continuous symmetry, 
low-temperature phases fall into the class of glass states which are known to be 
a difficult subject for theory. A usual theoretical approach in condensed matter 
physics is based on an exact solution of a simple related problem. Then one can 
develop e.g. a perturbative expansion to take into account the ingredients missed in 
the exactly solvable problem. A standard source of simple and useful exact solutions 
is mean field models. However, for strongly disordered systems the mean field 
approximation is often neither simple nor very useful. An example is the classical 
solution of the mean-field spin glass by Parisi.u If this very nontrivial result does 
capture physics of real short range spin glasses is still an open question. 

Much effort was devoted to theoretical understanding of strongly disordered 
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magnets, liquid crystals, vortex states of impure superconductors, and other related 
problems, but the lack of useful exactly solvable models has limited progress in 
the field. On the other hand, one could hope that for weakly disordered systems 
an appropriate solvable model is just an analogous system without disorder, and 
impurities can be taken into account with a perturbation theory. Apparently, the 
ability of weak randomness to destroy long range order, which is present in pure 
systems, is an argument against such an optimistic view. On a more rigorous level, 
the exact solution of the spherical model of amorphous magnets! - ! suggests that 
there is no qualitative difference between weakly and strongly disordered systems 
with continuous symmetry. Only recently has it been realized that this insensitivity 
to the disorder strength is an artifact of the spherical approximation, and a similarity 
between ordering in weakly disordered continuous systems and their pure analogs 
does exist. Long range order is prohibited in the systems of continuous symmetry 
in the presence of impurities, but instead quasi-long range order can emerge. This 
means that the average value of the order parameter over the volume is zero, but 
the correlation length is infinite and correlation functions obey power dependencies 
on the distance. In other words, long range order is only weakly broken. 

The concept of quasi-long range order (QLRO) was first introduced in the pure 
two-dimensional XY model.a Due to Abelian symmetry the degeneracy space of the 
XY model (i.e. the manifold each point of which represents a ground state) is locally 
indistinguishable from the degeneracy space of the free field. Thus, neglecting 
topological defects (this is legitimate at low temperatures) one can map the pure 
XY model onto the free field. This mapping gives an easy way to understand the 
low-temperature behavior of the XY model. The result is that long range order 
is present in spatial dimensions D > 2 and absent in lower dimensions. In two 
dimensions there is an intermediate situation of QLRO in the low-temperature 
phase. Our experience with pure systems shows that QLRO is a rare phenomenon. 
It is possible only in Abelian systems and only in the lower critical dimension that 
separates the regime of long range order in all higher dimensions from the regime 
with neither long nor quasi-long range order in all lower dimensions. 

Recent theoretical studies of impure systems0-|ll show however that QLRO 
can be much more common in the problems with disorder than in homogeneous 
systems. First, in a given system it can occur in a whole range of dimensions 
and not only at thelower critical dimension. Second, it is not prohibited by non- 
Abelian symmetry.L-3 E3 The possibility of QLRO was predicted in the Abrikosov 
state of impure superconductors,™ in amorphous magnetsE-Jiij and in uniaxial ne- 
matics confined in random porous matricesEj Although more work is needed to 
completely understand glass states of weakly disordered systems, and in particular 
numerical and experimental results are contradictory, we believe that this is only a 
small portion of the list of disordered systems which could possesa QLRO phases. 
Note that QLRO is possible only in weakly disordered systemsO The nature of 
glass states at strong disorder and disorder driven phase transitions from the weak- 
disorder to strong-disorder regime is an important open problem which is beyond 
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the scope of the present review. A very recent discussion o f th is question in the 
context of disordered superconductorstM can be found in Ref.La Note also that the 
existence of a QLRO state is sensitive to details of the systerrkFor example, it ex- 
ists in the random-anisotropy Heisenberg model but is absenirJ in the very similar 
random-field Heisenberg model that describes relaxor ferroelectrics.E3 

From the technical point of view the theoretical prediction of QLRO is based 
on the renormalization group applied to the weak-disorder fixed point. It turns out 
that an infinite set of relevant operators emerges in the problem and dealing with 
them requires some care. This technical problem takes its origin in the complicated 
structure of the energy landscape that contains a plethora of energy minima. The 
latter is certainly not surprising for a glass state. Earlv_attempts to apply the renor- 
malization group to this class of problems (e.g. Ref.E2h_failed due to an incorrect 
treatment of the complicated energy landscape: TheoryBj incorrectly predicted the 
dimensional reduction by 2 in comparison with corresponding pure systems and 
suggested that QLRO is possible only in 4-dimensional random systems of Abelian 
symmetry. 

In this paper we discuss recent theoretical results on the nature of glass states 
of weakly disordered systems with continuous symmetry. Since the problem of 
the Bragg glass state of the vortex lattice in impure superconductors was recently 
reviewed in excellent papers we consider impure superconductors only briefly 
and concentrate on systems of non-Abelian symmetry: disordered magnets and 
nematic liquid crystals. The central point is that in many cases the Larkin-Imry- 
Ma pictureHra of uncorrelated domains breaks down and quasi-long range order 
emerges. 

The article is organized as follows. In the next section we introduce models of 
some relevant systems: amorphous magnets, relaxor ferroelectrics, vortex lattices 
in disordered superconductors and uniaxial nematics in random porous media, and 
map them onto the random-field (RF) and random-anisotropy (RA) O(N) models. 
Section 3 contains a qualitative discussion of the glass phases. In section 4 we 
derive the functional renormalization group equations. Their solution is given in 
section 5. First we discuss a simple approximate solution and then proceed with 
a systematic approach. Some details of the systematic solution are summarized 
in the Appendix. In section 6 we compare the predictions with experiments and 
numerical experiments. In that section we also describe the rich phase diagram of 
disordered nematic. This phase diagram includes, in particular, two QLRO glass 
states. Section 7 contains conclusions. 

2. Disordered Systems of Continuous Symmetry 

In this section we introduce several disordered systems of continuous symmetry. 
Albeit physically very different they all can be described by the random-field or 
random-anisotropy O(N) models (or slight modifications of these models). Note 
that the problems of nematics and vortex lattices in superconductors can be mapped 
onto the 0(N) models in the low-temperature phases only. Mapping becomes in- 
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valid near the phase transitions to the high-temperature phases. 

2.1. Amorphous magnet 

If a solid is obtained with a fast freezing from a liquid the atoms may have not 
enough time to form a crystal lattice. In this case an amorphous solid is formed. 
One can imagine it as a liquid in which the positions of all particles are suddenly 
fixed. Similar to crystals amorphous solids have special directions in each point but 
these directions are different in different places and uncorrelated for distant points. 
In a magnetic system, special directions lead to the appearance of easy magnetic 
axes or planes. In contrast to crystals, in amorphous magnets the directions of 
easy axes are random. The simplest model aLamorphous ferromagnets is hence the 
random-anisotropy (RA) Heisenberg model.Ell It has the following Hamiltonian: 

H = -J^n l n,-Z?^(n,h l ) 2 , (1) 

(ij) « 

where J is the exchange strength, D the anisotropy strength, n; the spin in site 
i, {ij) denotes a pair of neighboring sites, and the unit vector describing the 
direction of random anisotropy in site i varies from site to site. 

At low temperatures amorphous magnets possess a spin-glass-like state. This 
fact can be easily understood as random anisotropy is strong since in this case the 
spins are frozen along their local anisotropy axis. A review on strongly disordered 
amorphous magnets can be found in Rcf.Ea In this paper we concentrate on the case 
of weak anisotropy. We shall see that in the latter case there are strong correlations 
between spins and the system possesses QLRO. An important consequence of this 
prediction is divergence of the magnetic susceptibility in the low-temperature phase. 

2.2. Nematic in aerogel 

Quenched disorder is inevitably present even in the most pure solids. This explains 
a lot of phenomena, e.g. the residual resistance of metals. On the other hand, liq- 
uids are usually homogeneous and introducing quenched disorder in them requires 
special efforts. One of the approaches consists in pouring a liquid into a randomly in- 
terconnected network of pores. Such liquid-porous-matrix systems emerge in many 
natural and technological processes giving rise to a lasting scientific activity. The 
recent surge oJLinterest in the field is due to a new micropore material: low-density 
silica aerogel.E^I Its density can be varied in a wide range up to more than 99% void 
volume fraction. This allows the investigation of both strongly and weakly con- 
fined fluids. The most interesting situation emerges in systems with many degrees 
of freedom, e.g. He-3 c3 and liquid crystals £3 _ Eil In these substances the porous 
matrix not only geometrically confines the liquid but also induces a random ori- 
enting field that fixes the direction of the order parameter near the surface of the 
matrix. The random-field disorder is known to cause spin-glass effectsEa and such 
phenomena were indeed observed experimentally in liquid-crystal-aerogel systems. 
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In particular, a slow glassy dynamics was reported in Refs.EE 

Our aim is to describe the glass state of nematics in random porous media in 
the case of weak disorder. This corresponds e.g. to a weak interaction with the 
surface of a porous matrix. We postpone the discussion of the disorder strength to 
section 6 and formulate the simplest model of disordered uniaxial nematics. This 
is a slight complification of the RA Heisenberg model (Q) . The free energy density 

F = F d + F pm (2) 
of the nematic in the porous matrix includes the Frank distortion energyil 

F d = [/^(divn) 2 + ^(ncurln) 2 + K 3 (n x curln) 2 ]/2, (3) 

where n is the director, and the interaction F pl[L with the surface of the random 
matrix. The interaction tends to align the director parallel to the surfaceO We 
model the interaction as 

F pm = (hn) 2 , (4) 

where h is a random vector representing the normal to the surface. This is the 
simplest choice compatible with the equivalence of the opposite orientations of the 
director. Due to the universality the model captures all large-scale physics. The 
average amplitude of the random vector h is a measure of the disorder strength. It 
is a phenomenological parameter which depends on the pore size, anchoring energy 
and fractal structure of the porous matrix. 

We demonstrate that the model ((§-[§) possesses QLRO in its low-temperature 
phase at weak disorder. As a consequence the light-scattering cross-section diverges 
at small angles. We also consider effects of external electric and magnetic fields, and 
mechanical deformations of the porous matrix. It turns out that at weak disorder 
there are several phases including two glass states. 

2.3. Relaxor ferroelectric 

Relaxor ferroelectrics are interesting because of their diffuse phase transitions ex- 
tending over a finite range of temperatures. Most relaxor ferroelectrics are dis- 
ordered ionic str uct ures, in particular, solid solutions. In the best known re- 
laxor ferroelectrkip PbMg 1/ / 3 Nb 2 /303 the spontaneous polarization vanishes at 
all temperatures^ and the phase transition is believed to be destroyed by the 
randomness. E3 The randomness can be described as quenched random electric fields.E3 
This leads to the random-field Heisenberg model of the relaxor ferroelectric. Its 
Hamiltonian reads 

H = - ^ JijTH-a-j - ^(ruhi), (5) 

i 

where is the local dipole moment and the random electric field. This model 
is certainly a simplification. In particular, it misses nonrandom anisotropy (which 
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is cubic in PhMg^N^/sOs). Nonrandom anisotropy leads to a discrete symmetry 
group and hence can stabilize long ra nge (ferroelectric) order which is destroyed by 



random fields. □ However, it is expectecO that in PbMg!/ 3 Nb2/303 the anisotropy is 
weak. We demonstrate that the model (^) with the interaction of the nearest neigh- 
bors only does not possess QLRO. The effect of the long range dipole interactions 
which are present in relaxor ferroelectrics is an interesting open question. 

2.4. Vortices in disordered superconductor 

Disordered superconductors were recently considered in excellent reviewsBE so in 
this subsection we only provide a basic background and briefly discuss mapping 
onto the random-field XY model. Our discussion of that model in the subsequent 
sections is also brief. More details can be found in the abovementioned reviews. 

In the absence of disorder, vortices in type II superconductors form a lattice. 
Small fluctuations of the vortices about their equilibrium positions can be described 
with elastic theory. Let us denote by Ri the equilibrium position of the vortex in 
lattice site i. The vortex displacement u(i?;, z) = u l (z) relative to its equilibrium 
position is a two-component vector. In the absence of dislocations we can assume 
that the displacement field is slowly varying on the scale of the lattice spacing 
a. We will see in section 5 that topological defects are indeed irrelevant at weak 
disorder and low temperatures. Hence, a continuous description can be used. The 
continuous elastic Hamiltonian has the following form 

1 v f d 3 q 

ff el=2 2j J J^3®M<l) U a(q)uf3(-q), (6) 

where a, j3 = x,y label the coordinates, and the structure of the matrix Q a p for 
the triangular Abrikosov lattice can be found in Ref.EH For simplicity we consider 
below a one-component field u. This corresponds to an anisotropic superconductor 
such that the vortices can be displaced in one direction only. It turns out that this 
simple model allows to obtain a good quantitative description of some aspects of 
the large-scale behavior of disordered isotropic superconductors. A discussion of a 
more accurate model taking into account two components of the displacement field 
and the triangular symmetry of the Abrikosov lattice can be found in RefsSS We 
briefly discuss the results of Refs@'@ in section 6. The following elastic Hamiltonian 
describes the one-component displacement field: 

H ^ = \J d 3 x(Vuf. (7) 

This is the Hamiltonian of the XY model without dislocations. 

The next step is to include disorder. The pinning energy readsO 

H pin = J d 3 xp(x.)V(x), (8) 

where the vortex density 
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p(x) =£«S(x-R i -u i ), (9) 

i 

and V is the random pinning potential with the following statistics: 

(V(x)) = 0; (V(x)V(y)} = A<5(x - y). (10) 

Mapping on the random field XY model is little bit tricky. We rewrite the vortex 
density in the following way 

p(x) = J2 J d 3 x'6(x - x' - u i (x'))5(R i - x') 

i 

= po / d 3 iz/£(x — x 7 — exp(zQx') 

Q 

= podet -1 ^^ + <9 ct u /5 ]^exp(iQ[x- u]) 

Q 

~ po[l - d a u a ] ^exp(iQ[x- u]), (11) 

Q 

where po is the average vortex density, Q are reciprocal lattice vectors. Substitut- 
ing ( |ll| ) into (|J) we get an infinite set of random contributions to the energy of the 
form po[l — 9q,m q ]Vq(x) exp(— iQu(x)), where Vq are random fields. We are inter- 
ested in the large distance behavior and hence keep only the most relevant terms 
which should be taken into account by the renormalization group. Thus, we neglect 
derivatives of u. We shall see that an infinite set of other random contributions 
is relevant in the RG sense. Here we only write the minimal Hamiltonian of the 
correct symmetry for a one-component field u: 

H = ^ J d 3 x(Vu) 2 + J d 3 xh(x) cos(tt(x) - a(x)), (12) 

where h and a are random. This is the Hamiltonian of the random-field XY model. 

It turns out that a QLRO state (Bragg glass) emerges in impure superconductors 
at weak disorder and low temperatures. The name 'Bragg glass' is due to sharp 
Bragg peaks which are observed in this state because long range order is only weakly 
broken. Note that the random field h depends on the vortex density po and hence on 
the applied magnetic field along which the vortices are directed. Since the density is 
inversely proportional to the magnetic field, the effective disorder strength is small 
at weak magnetic fields. Indeed, in a given sample (i.e. at a fixed pinning potential 
V(x)) the Bragg glass state is observed at low magnetic fields. 

3. Glass States 
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In this section we qualitatively discuss glass states of systems with continuous 
symmetry. We demonstrate that disorder, which breaks the continuous symme- 
try, makes long range order impossible and discuss a possibility of QLRO. A more 
general question is the existence of a phase transition to a glass state with or without 
QLRO. There is a simple argument in favor of such transition in weakly disordered 
systems whose pure analogs have a first order transition between high- and low- 
temperature phasesO Note that in homogeneous nematics the phase transition is 
of the first order. The phase transition from vortex liquid to vortex solid in pure 
superconductors is also believed to have first order .0 

First, let us show that even at arbitrarily weak disorder there is no long range 
order in disordered systems of continuous symmetry. We employ the Imry-Ma 
argument .□ It can be used in any system in which disorder breaks the continuous 
symmetry, but we restrict our discussion by the random-field and random- anisotropy 
O(N) models which are considered in more details in the following sections. 

We start from the the classical nonlinear er-model with the Hamiltonian 

H= /d D x[J^^n(x)^n(x) + V5 mp (x)], (13) 
J n 

where n(x) is the unit vector of the TV-component magnetization, Vi mp (x) the ran- 
dom potential. In the RF case it has the form 

V imp = -^Mx)n a (x); a = l,...,N, (14) 

a 

where the random field h(x) has a Gaussian distribution and {h a (x)hp(x')) = 
A 2 S(x — x')6 a p. In the RA case the random potential is given by the equation 

Vimp = -y^r a;3 (x)n a (x)ng(x); a, (3 = 1, N, (15) 

where r Q ^(x) is a Gaussian random variable, (r QI| g(x)r 7< j(x / )) = A 2 8 ai 5ps5{'x. — 
x'). The random potential Jl5| ) corresponds to the same symmetry as the more 
conventional choice Vi mp = — (hn) 2 used in section 2 but is more convenient for the 
further discussion. 

We assume that the temperature is low and thermal fluctuations are negligible. 
The Imry-Ma argumentBcil suggests that in our problem long-range order is absent 
at any dimension D < 4. Let us assume that the magnetization changes by order 
1 at scale L. This costs the exchange energy of order E e ~ JV/L 2 , where V is 
the volume. It should be compared with a possible energy gain due to disorder 
potential E d ~ AV/L D x L D / 2 = AV/L D ' 2 . One finds that at D < 4 the exchange 
energy E e < Ed as L is large. Hence, it is favorable to destroy long range order. 

One can estimate the Larkin length, up to which there are strong ferromagnetic 
correlations, with the following qualitative RG approach. Let one remove fast modes 
and rewrite the Hamiltonian in terms of the block spins, corresponding to the scale 
L = ba, where a is the ultraviolet cut-off, b > 1. Then let one make rescaling such 
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that the Hamiltonian would restore its initial form with new constants A(L), J{L). 
Dimensional analysis provides estimations 



J(L) ~ b D - 2 J(a); A(L) ~ b D/2 A(a). (16) 

To estimate the typical angle <f> between neighboring block spins, one notes that the 
effective field, acting on each spin, has two contributions: the exchange contribution 
and the random one. The exchange contribution of order J(L) is oriented along the 
local average direction of the magnetization. The random contribution of order A(L) 
may have any direction. This allows one to write at low temperatures that </>(£) ~ 
A(V)I J(L). The Larkin length corresponds to the condition <fr(L) ~ 1 and equals 
L ~ (J I A) 2 I ( 4 ~-°) in agreement with the Imry-Ma argument.0 If Eq. ( [l6| ) were exact 
the Larkin length could be interpreted as the correlation length. However, there 
are two sources of corrections to Eq. (|l^) . Both of them are relevant already at the 
derivation of the RG equation for the pure system in 2 + e dimensions .e3cj The first 
source is the renormalization due to the interaction and the second one results from 
the magnetization rescaling which is necessary to ensure the fixed length condition 
n 2 = 1. The leading corrections to Eq. ( |l6| ) are proportional to tfi 2 J,(j) 2 A. Thus, 
the RG equation for the combination (A(L)/J(L)) 2 reads 

dlnL\j(L)J e {j(L)J +C \J(L)J ' [U) 

If the constant c in Eq. ([l7|) is positive the Larkin length is the correlation length 
indeed. But if c < the RG equation has a fixed point, corresponding to the phase 
with an infinite correlation length. As seen below, both situations are possible, 
depending on the number N of the magnetization components. The case of the 
infinite correlation length corresponds to QLRO. We shall see that QLRO is possible 
as N < N c , where iV c is a critical number. This critical number is larger in the 
RA model, since the fluctuations of the magnetization are stronger in the RF case. 
Indeed, in the RF model the magnetization tends to be oriented along the random 
field, whereas in the RA case there are two preferable local magnetization directions 
so that the spins tend to lie in the same semispace. 

The above scaling argument does not allow to understand if the low-temperature 
phase is different from the high-temperature paramagnetic state in case when QLRO 
is absent. The answer turns out to be positive, at least if disorder is weak and the 
system would have a first order phase transition in the absence of disorder .c3 

Near a first order phase transition a pure system has two phases: a disordered 
phase with a finite correlation length and an ordered one with an infinite correlation 
length. In the presence of impurities the ordered state is broken into a set of Imry- 
Ma domains of size L ~ {J/A) 2 / {A - D \ This lowers the free energy that would equal 

Fi(T) in the pure system by E x A(V/L D )L D / 2 VA 4 /^ D \ As disorder is 

weak the correlation length I of the disordered phase is less than the Larkin length L 
and we can estimate the free energy gain in the presence of disorder as E% ~ —VxA 2 , 
where \ is the susceptibility. The condition of the first order transition is that the 
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free energies of two phases Fx (T) + E\ and F2 (T) + E2 are equal. One finds the 
following shift of the phase transition temperature 



^ ■jw-fV (18) 

In 3 dimensions we obtain AT ~ A 2 . This argument does not work for systems 
with second order transitions since near the phase transition the correlation lengths 
of both phases are infinite. It is also valid only at D > 2. In two dimensions first 
order transitions are prohibited in disordered systems: They are destroyed due to 
roughening of domain walls : lj 

The last remark is that QLRO is not expected in strongly disordered syste 
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This agrees with the structure of the phase diagram of impure superconductors 
For the RF O(N) model the absence of QLRO at strong disorder is proven rigorously.llll 

4. Functional Renormalization Group 

In the previous section the RG equations are discussed from the qualitative point of 
view. Eq. ( |i"7|) corresponds to the simple Migdal-Kadanoff approach of the section 
5.1. In the present section we derive the RG equations in a systematic way. 

The one-loop RG equations fox the iV-component RF and RA models in 4 + e 
dimensions were derived in Ref.E3 Their solution allowed to describe the phase 
transition from the .ferromagnetic to paramagnetic state of the RF O(N) model 
above 4 dimensions .c3 The RG equations in dimensions D < 4 can be obtained by 
just changing the sign of e in the RG equations in_D > 4. 

In our derivation of RG equations we follow Ref.li3 We use the method, suggested 
by Polyakovocj for the pure system in 2 + e dimensions. This method is technically 
simpler and closer to the intuition than the other approaches. A disadvantage of 
the method is the difficulty of the generalization for_the higher orders in e. This 
generalization requires the field-theoretical approach.cJ 1—11—1 

The same consideration as in the and random manifoldE§E2 models sug- 

gests that near a zero-temperature fixed point in 4 — e dimensions Ihcre is an infinite 
set of relevant operators. Let us show that after replica averaging^ the relevant part 
of the effective replica Hamiltonian can be represented in the form 



H, = I d D 



.r 



V — V9 n d n - V R ^ anb) 

. a /j ab 



(19) 



where a, b are replica indices, R(z) is some function, T the temperature. We ascribe 
to the field n the scaling dimension 0. We also assume that the coefficient before 
the gradient term in ( |l9| ) is 1/(2T) at any scale. Then in the (4 — e)-dimensional 
space the scaling dimension of the temperature At = —2 + O(e). Any operator 
A containing m different replica indices is proportionalcj to 1/T m . Indeed, before 
replica averaging any term of the Hamiltonian contains one replica index and the 
temperature in the minus first power. If we expand the partition function in Taylor 
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series any m-replica term of the expansion contains T~ m . This property conserves 
after disorder averaging. To obtain the effective replica Hamiltonian one reexpo- 
nentiates the disorder-average series. Then one easily sees that the power of 1/T 
in anv term of the replica Hamiltonian equals the number of the replica indices 
in it.Ell Hence, the scaling dimension A^ of the operator A satisfies the relation 
A^ = 4 — n + mAx + 0(e), where n is the number of the spatial derivatives in 
the operator. The relevant operators have A^ > 0. Hence, the relevant operators 
cannot contain more than two different replica indices. A symmetry consideration 
shows that all possible relevant operators are included into Eq. (|l9|). The function 
R(z) is arbitrary in the RF case. In the RA case R(z) is even due to the symmetry 
with respect to changing the sign of the magnetization. 

The one- loop RG equations for the TV-component model in 4 —fj, dimension S_arfi 
obtained by a straightforward combination of the methods of Ref.E3 and Refs.e3e3 
We express each replica n°(x) of the magnetization as a combination of fast fields 
0°(x),i = 1, N — 1 and a slow field n' a (x) of the unit length. We use the 
representation 



n a (x)=n' a (x) l-^(x))H^^(xK(x), (20) 

y i i 

where the unit vectors e°(x) are perpendicular to each other and the vector n' a (x). 
The slow field n' a can be chosen in different ways. For example, one can cut the 
system into blocks of size L>a, where a is the ultra-violet cut-off. In the center 
x of a block the vector n' a (x) should be parallel to the total magnetization of the 
block. In other points the field n' a should be interpolated. We assume that the 
fluctuations of the magnetization are weak, that is (<j>f) <C 1. Then the fluctuations 
of the field n a are orthogonal to the vector n' a because of the fixed length constraint 
(n a ) 2 = 1. 

To substitute the expression (|2(]) into the Hamiltonian we have to differentiate 
the vectors e". We use the following definition 

!^ = -c>'«+£/^. (21) 

" k 

It is easy to showSIl that ^i( c %i) 2 — X^(^ n ' a ) 2 - With the accuracy up to the 
second order in (j) the replica Hamiltonian dl^) can be represented as follows 



Ho = I d 



D , 



^ J2{R(n' a n' b ) + A ab (tff + Bgfifl + C##$}j , (22) 



ab 



where the summation over the repeated indices i,j,k,fi is assumed and 
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A ab = _( n 'a n "')i ? '( n '° n "») ; ^ = ( n ' b ef)(n ,b e])R"(n' a n' b ); 

C? b = {e1e b )R'{n' a n' b ) + {n' a e b )(n' b e?)R" (n'V 6 ). (23) 

In the last formula R' and R" denote the first and second derivatives of the function 



R(z). We have omitted the terms of the first order in cj> in Eq. (22). These terms 
are proportional to the products of the fast field cf> and some slow fields. Hence, 
they are non-zero only in narrow shells of the momentum space. One can show that 
their contributions to the RG equations are negligible. 

To obtain the RG equations we have to integrate out the fast variables (j>f . Near 
a zero-temperature fixed point the Jacobian of the transformation n — > (n', 4>i) can 
be ignored, since the Jacobian does not depend on the temperature. The integration 
measure is determined from the condition that the field n' a is a slow part of the 
magnetization. This condition imposes restrictions on the fields (f>. The expression 
( p2| ) depends on the choice of the vectors e° (|20|). However, after integrating out the 
fields (j> the Hamiltonian can depend only on the slow part n' a of the magnetization. 
One can make the calculations simpler, considering special realizations of the field 
n' a . To find the renormalization of the disorder-induced term R(z) ( |l9| ) we can 
assume that the field n /a does not depend on the coordinates. The renormalization 
of the gradient energy can be determined, assuming that the vectors n' a (x) depend 
on one spatial coordinate only and lie in the same plane. In both cases the vectors 
e° can be chosen such that in Eq. ( pl| ) /" ik = 0. At such a choice the integration 
measure can be omitted and the fields (f>" can be considered as weakly interacting 
fields with the wave vectors from the interval 1/a > q > 1/L. 

To derive the one-loop RG equations we express the free energy via the Hamilto- 
nian (^2|) . Then we expand the exponent in the partition function up to the second 
order in {Hr — J d D x Y^ni^^i) 2 / i^T)) and integrate over <f>f. Finally, we make 
a rescaling. The vectors e" can be excluded from the final expressions with the 
relation X)i( ae ?)(be") = (ab) — (an' a )(bn' a ). In a zero-temperature fixed point 
the RG equations are 



dlnT 



-{D-2) + 2(N -2)R'(l) + 0{R 2 ,T); (24) 



d\nL 

= = eR(z) + 4(N - 2)R(z)R'(l) - 2{N - l)zR'(l)R'(z) 

a In L 

+2(1 - z 2 )R'(l)R"(z) + (R'(z)) 2 (N ~2 + z 2 ) 
-2R'(z)R"(z)z(l - z 2 ) + {R"{z)) 2 {\ - z 2 ) 2 , (25) 

where the factor l/(87r 2 ) is absorbed into R(z) to simplify notations. The RG 
equations become simpler after the substitution of the argument of the function 
R(z): z = cos0. In terms of this new variable one has to find even periodic 
solutions R(4>). The period is 2ir in the RF case and tt in the RA case due to the 
additional symmetry of the RA model. The one-loop equations get the form 
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dhxT 
d\nL 



(D-2)-2{N-2)R"(0) + O(R 2 ,T); 



(26) 



u = = eR(<f>) + (R"m 2 - 2R"(<f>)R"(0)- 



(N-2) 



d\nL 



4R(c/>)i?"(0) + 2ctg(f)R , (<j))R"(0) 



R'{<t>) 



0(R 3 ,T) (27) 



Eq. ( Pq ) provides the following result for the scaling dimension At of the temper- 
ature 



-2 + e - 2(N - 2)R"(0). 



(28) 



The two-spin correlation function is given in the one-loop ordei@'@ by the ex- 
pression 



<n"(xK(x')) = (n'»(x)n"V)) (l - (^W) 3 



Hence, in the fixed point (n(x)n(x')) ~ |x — x'| n , where 



rj = -2(N -l)R"(<f> = 0) 



(29) 



(30) 



Let us find the magnetic susceptibility in the weak uniform external field H . We 
add to the Hamiltonian © the term - £a / d D xHn a z /T (the field is directed along 
the z-axis) . The renormalization of the field H is determined by the renormalization 
of the temperature ( ^rj ) and the field n. In the zero-loop order the renormalized 
magnetic field h(L) depends on the scale as h(L) = H x (L/a) 2 . Hence, the corre- 
lation length R c ~ H^ 1 / 2 . The magnetization, averaged over a block of size R c , is 
oriented along the field. The absolute value of this average magnetization is propor- 



tional to R, 



-n/2 



This allows us to calculate the critical exponent 7 of the magnetic 



susceptibility x{H) ~ H 7 in a zero-tempcrature fixed point: 



■y = l + (N-l)R"(<f) = 0)/2. 



(31) 



In Rcf.EJ Eqs. ( J24| , |25| ) were derived with a different method. In that paper the 
critical behavior in 4 + e dimensions was studied by considering analytical fixed 
point solutions R(z). In the Hejsenberg model, analytical solutions are absent and 
they are unphysical for N ^ 3H3 In 4— e dimensions appropriate analytical solutions 
are absent for any TV. To demonstrate this let us differentiate Eq. (^5|) over z at 
z = 1. For any analytical R(z) we obtain the following flow equation 



dR'jz = 1) 
dlnL 



eR'{z= 1) + 2(N -2){R'(z = l)) 2 



(32) 
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At TV > 2 the fixed point of this equation R'(z = 1) = -e/[2(N - 2)] < 0. It 
corresponds to the negative critical exponent rj ( p0| ) and hence is unphysical (two- 
spin correlation functions are limited and cannot grow up to infinity as R — > oo). 
However, we shall see that in the RA model some appropriate non-analytical fixed 
points R(z) appear. In these fixed points R"{z = 1) = oo. In Rcf.E3 the RG charges 
are the derivatives of the function R(z) at z = 1. Thus, in a non-analytical fixed 
point these charges diverge. In the systems with a finite number of the charges their 
divergence implies the absence of a fixed point. However, if the number of the RG 
charges is infinite such a criterion does not work and is even ambiguous. Indeed, 
the set of charges can be chosen in different ways and e.g. the coefficients of the 
Taylor expansion about z = remain finite in our problem. 

Nonanaliticity of the fixed point solution is related with a complicated structure 



of the energy landscape 



nxcc 

BE 



interested reader may refer to Ref.0 



This has important consequences for dynamics. An 



5. Solution of the Renormalization Group Equations 

In this section we solve the RG equations of the preceding section. The main aim 
is to understand in which cases there is QLRO and when it is absent. The RG 
equations are the same for the RF and RA O(N) models. The difference between 
the models is the symmetry of the fixed point solution. In the RF case we are 
looking for the solution which is stable to arbitrary perturbations and describes 
the low-temperature phase of a weakly disordered system. In the RA case one has 
to respect the symmetry to changing the sign of the magnetization. As discussed 
above in the RA case the function R{z) must be even. The solution should be stable 
only to the perturbations that do nut break this property. 

We obtain the following resultaLa : In the RF model, QLRO exists at N = 2 
only. In the RA case, QLRO is present in 4 — e dimensions for N < 10. The critical 
exponents describing the QLRO phases of the RA O(N) model are given in Table 
1. 

The section is organized as follows: In the first subsection we discuss a simple 
approximate solution of the RG equations. In two next subsections we provide a 
systematic analysis of the RF and RA models. Since our RG approach ignores 
topological defects we have to check that they are irrelevant. This is done in the 
fourth subsection. 



5.1. Migdal-Kadanoff renormalization group 

This subsection contains a simple approximate version of the renormalization group. 
The results for the critical exponents of the XY and Heisenberg models have a very 
good accuracy. The value of the magnetization component number N c , at which 
QLRO disappears in the RF model, is probably exact. However, the critical number 
of the components in the RA model is underestimated. 



5.1.1. Random field 
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We use the following ansatz for the disorder-induced term in the Hamiltonian (|l9|): 
i?(n a n{,) = an a rib + /3, where a and (3 are constants. This expression corresponds 
to the Gaussian RF disorder (|l4|). Below we ignore the generation of other con- 
tributions to the function R(z). The missed contributions are related to random 
anisotropics of different orders. In terms of the angle variable i 



R(<f>) = acos<t> + p. (33) 



To ensure consistency we have to truncate the RG equation (|27|) . We substitute 
the ansatz (|33|) into Eq. ( p7| ) and retain only the terms, proportional to cos0 or 
independent of 0. This leads to the following RG equation for the constant a ( |33| ) 

= ea + 2a 2 (N-3). (34) 

cl In Li 

For N < 3 Eq. (|34| ) has a stable solution a = e/[2(3— AT)]. This solution corresponds 
to a QLRO state. The critical exponent (|30| ) equals 

_ (N-l)e 

(3 — TV) (d5j 

At AT = 2 this result has less than ten percent difference from the systematic 
theory.ta QLRO disappears at A^ = 3. This is the same critical number which is 
found with a systematic approach. 

For A^ > 3 a fixed point exists in 4 + e dimensions. It describes the transition 

he critical 
3 However, 



between the ferromagnetic and paramagnetic phases. In this fixed poin 
exponent ( |35| ) satisfies the modified dimensional reduction hypothesis, 
this is an artifact of ihe Migdal-Kadanoff approximation, since the correct value of 
the critical exponento differs form Eq. ( |35| ) . The detailed discussion of the critical 
exponents of the RF O(N) model in 4 + e dimensions is beyond the scope of the 
present article. Some details can be found in Ref.c3 

5.1.2. Random anisotropy 

In this case we use the ansatz R(n a THb) = A(n a rib) 2 + B. In terms of the variable 
(j) Eqs. ( p6j^7| ) R(4>) — a cos 2(f) + [3. We again substitute our ansatz into Eq. ( |27| ) 
and retain the terms, proportional to cos 20, and the terms, independent of <fi. The 
RG equation for the constant a has the form 

<!n = ea + 8(N-6)a 2 . (36) 



dlnL 

The fixed point solution of this equation is a = e/[8(6 — N)]. It describes the QLRO 
phase at N < 6. At N = 3 the function R((f>) = a cos 20 + (3 is just i? w= o of section 
5.3.3. The critical exponent of the two-spin correlation function is given by the 
following equation 

V = T^AT- (37) 
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At N — 2, 3 this value is close to the results of the systematic approach (Table 1). 

The approximate analysis suggests that QLRO disappears at N = 6. This is 
different from the result of the systematic approach N c = 10. 

5.2. Random field 

5.2.1. N = 2 

For the RF XY model the one-loop RG equations ( p6| , ^7|) can be solved analytically.El 
The solution is a periodic function R(4>) with period 2n. In interval < 4> < 2ir the 
fixed point solution R(cf>) is given by the formula 



?r 4 e 



R{<t>) = Y L 1 / 36 ~ W*/ 2 ^ 2 ~ (^)) 2 J ■ ( 38 ) 

This is a stable fixed point. This can be verified with the linearization of the flow 
equation ( |27j ) for small deviations from the fixed point. 

The solution j3q ) corresponds to QLRO with the critical exponents ( |30| , |3l| ) 
i] = 7r 2 /9e,7 = 1 — 7r 2 /18e. In the first order in e the exnonent i] equals the 
prefactor C before the logarithm in the correlation functionE 2 ] of the angles <^>(x) 
between the spins n(x) and some fixed direction: ((0(xi) — </>(x2)) 2 ) = Cln|xi— X2|. 
We expect that this coincidence does not extend to the higher orders. 

5.2.2. N > 2 

If N 7^ 2 the RG equation ( |27| ) is more complex. Fortunately, at N > 3 there is still 
a simple method to study Jiie large-distance behavior. The method is based on the 
Schwartz-Saffer inequality^ and shows that QLRO is absent. 

In Ref£3 the inequality is proven for the Gaussian distribution of the random 
field. It can also be proved for the arbitrary RF distributionJla We prove the 
inequality in the Appendix. 

Let us demonstrate the absence of physically acceptable fixed points in the RF 
case at N > 3. We derive some inequality for critical exponents. Then we show 
that the inequality has no solutions. Wfiuse a rigorous inequality for the connected 
and disconnected correlation functions^ 



((n(q)n(-q))) = <n Q (q)n Q (-q)) 

- (n a (q)n fc (-q)) < const y/ (n a (q)n a (-q)), (39) 

where n(q) is a Fourier-component of the magnetization, a, b are replica indices. The 

disconnected correlation function is described by the critical exponent (|3C|). The 
large-distance behavior of the connected correlation function in a zero-temperature 
fixed point can be derived from the expression 

X ~ /«n(0)n(x)))^ (40) 
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and the critical exponent of the susceptibility (|3l]). The integral in the right hand 
side of Eq. ([fOj) is proportional to R®~ m , where R c is the correlation length in the 
external field H, rji the critical exponent of the connected correlation function. For 
the calculation of the exponent 7 (|3l]) we used the zero-loop expression of R c via 
H. Now we need the one-loop accuracy. In this order R c ~ jy-VP-C-W-^-R 
This allows us to get the following equation for the exponent rji 

r]i=D-2- 2R"(0). (41) 

In a fixed point Eq. ( |39|) provides an inequality for the critical exponents of the 
connected and disconnected correlation functionsS The inequality has the form 

2(2- D + rjt) >i-D + n. (42) 
This allows us to obtain the following relation 

3 - iV 

4-£»< ? - TJ7 + (E), (43) 

where 77 is given by Eq. ( J30| ) . The two-spin correlation function cannot increase up 
to infinity as the distance increases. Hence, the critical exponent 77 is positive. At 
N > 3 this is incompatible with Eq. (Q) at small e. Thus, there are no appropriate 
fixed points for N > 3. This suggests the strong coupling regime with a presumably 
finite correlation length. We see that QLRO disappears at N > N c < 3. Numerical 



analysis of the RG equations supp 
Certainly, in the RF XY mod£ 
fixed points of the RG equation 



N c <3. 

Eq. ( f43| ) is satisfied. However, the unstable 
do not satisfy the inequality. 



5.3. Random anisotropy 

In this subsection we investigate a possibility of QLRO in the RA O(N) model. The 
first subsubsection is devoted to the simplest case of the XY model. The second 
subsubsection contains an inequality for the critical exponent rj. The derivation of 
the inequality is analogous to the derivation of Eq. (fl3|). This inequality is applied 
in the next subsubsections. The third subsubsection contains the results for the 
Heisenberg model. In the last subsubsection we consider the case N > 3. 

5.3.1. N = 2 

This case is studied analogously to the RF XY modellll At N — 2 the RG equation 
( p?] ) can be solved analytically. Its solution is a periodic function with period it. In 
interval < < 7r the fixed point solution R(<f>) is given by the formula 

m = m [ 1/36 _ (0A)2 (1 _ W7:))2 ] ■ (44) 

It is a stable fixed point. This can be verified with the linearization of the flow 
equation (p^) for small deviations from the fixed point. Another proof of stability 
is based on the inequality of the next subsection. 
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The stable fixed point corresponds to the QLRO phase at low temperatures and 
weak disorder. The critical exponents rj — 7r 2 e/36,7 = 1 — 7r 2 e/72. 

The solution ( f44| ) is non-analytical at (f> = 0, since R IV (cf> = 0) = oo. Hence, the 
Taylor expansion over <j) is absent. However, a power expansion over \cj>\ exists. We 
shall see below that the same behavior at small <j> conserves also at other N. 

5.3.2. An inequality for a critical exponent 

We use the same approach as in the RF model. Since in the RA case the random field 
is conjugated with a second order polynomial of the magnetization, the Schwartz- 
Soffer inequality!!] should be applied to correlation functions of the field m(x) = 
(n z (x)) 2 — 1/N, where n z denotes one of the magnetization components, 1/N is 
subtracted to ensure the relation (m) = 0. 

To calculate the critical exponent \x of the disconnected correlation function we 
use the representation ( |2p| ) and obtain the relation 

(™ a (x)™ a (x')) = (m' a (x)m'°(x')) ^1 - ^^^f ^ ) , (45) 

where a is a replica index, m' — (n' z ) 2 — 1/N the slow part of the field m. One finds 
H = -ANR"(0). 

The critical exponent /zi of the connected correlation function is determined 
analogously to the RF case. We apply a weak uniform field H, conjugated with 
the field m, and calculate the susceptibility dm/dH in two ways. The result for the 
critical exponent is /ii = D - 2 - 2(N + 2)R"(0). 

The Schwartz-Soffer inequality provides a relation between the exponents fj, and 
fii . It has the same structure as Eq. ( ft2] ) . Finally, we obtain the following equation 

V>^^-(N-l) + o(R). (46) 

In terms of the RG charge R(<j>) this inequality can be rewritten in the form 

R"(0) < -e/8 + o(R). (47) 

As discussed in Ref.0 any solution of the RG equations that does not satisfy 
Eq. ([IT]) describes an unstable fixed point. 

5.3.3. N = 3 

In this case we solve Eq. ( p7| ) numerically. Since coefficients of Eq. (|27|) are large as 
<^> — *■ 0, it is convenient to use a series expansion of the fixed-point solution i?(</>) at 
small </>. At larger <fi the equation can be integrated with the Runge-Kutta method. 
The following expansion over t — yj (1 — z)/2 = \ sin(</>/2)| holds 



____ +2asln -±— y — |sm -| 
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2a 



sin 4 Z + OQ sin 5 ^|), (48) 



3 3(N + A)J 2 vl 2 

where a = R"(<fi = 0)/e. We see thajLthe RG charge R(<f>) ismpn-analytical at small 
</>. Similar to the random manifoldEiEJ and random-field XYI13 models R IV (0) = 



oo. 



Numerical calculations show that at any TV the solutions, compatible with the 



inequality fl47j), have sign "+" before the third term of Eq. (48). The solutions 



to be found are even periodical functions with period tt. Hence, their derivative is 



zero at tj> = tt/2. At N — 3 there is only one solution that satisfies Eq. (47). It 
corresponds to 

R"(4> = 0) = -0.1543e. (49) 

If this solution is stable Eqs. ( |30| , pl| ) provide the following results for the critical 
exponents 

77 = 0.626; 7= 1 - 0.15e. (50) 



All the other solutions of Eq. fl27| ) do not satisfy Eq. ( |47| ) and hence are unstable. 

We have still to test the stability of the solution found. For this aim we use an 
approximate method. First, we find an approximate analytical solution of Eq. (|27j). 
We rewrite Eq. @, substituting u(R"(4>)) 2 for (R"(4>)) 2 . The case of interest is 
u> = 1 but at lj = the equation can be solved exactly. The solution at w = 1 can 
then be found with the perturbation theory over u>. The exact solution at u> = 
is Ruj=o{4>) = e(cos 20/24 + 1/120). The corrections of order io k are trigonometric 
polynomials of order 2(k + 1). The first correction is 

Ri{4>) = — ^tt cos20 + cos40 + const (51) 
99 264 

After the calculation of the corrections we can write an asymptotic series for the 
critical exponent r\ (|30|): r\ = e(0.67 — 0.08cj + 0.14o; 2 — . . .). The resulting estimation 
r\ = e(0.67±0.08) agrees with the numerical result ( |50| ) well. This allows us to expect 
that the stability analysis of the solution i? w= o of the equation with lo = provides 
information about stability of the solution of Eq. (p7|). 

To study stability of the exact solution of the equation with ui = is a simple 
problem. We introduce a small deviation r (</>): R{<t>) = Ruj=q{4>) + r((f>) and write 
the flow equation for this deviation: 

= (5rU) + r"(<p) + r"(0) cos20)/3 + const x r"(0). (52) 

am L 

It is convenient to use the Fourier expansion r(0) = ^2 m a m cos 2m<j). The flow 
equations for the Fourier harmonics can be easily integrated. We see that a m — > 
as L — > oo for any m > 0. The solution is unstable with respect to the constant 
shift dQ. However, this instability has no interest for us, since correlation functions 
do not change at such shifts.c2l Indeed, a constant shift corresponds to the addition 
of just a random term, independent of the magnetization, to the Hamiltonian (|l3j). 
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Thus, the RG equation possesses a stable fixed point. This fixed point describes 
the QLRO phase with the critical exponents (p0[). 

As usual in critical phenomena, in 4 dimensions the one-loop RG equations allow 
one to obtain the exact large-distance asymptotics of the correlation function. In 
the 4-dimensional case R{4>) = &((/>)/ In L, where R{4>) satisfies Eq. ( gTp at e = 1. 
We obtain the following result for the two-spin correlation function with Eq. ( p9| ) 

(n(x)n(x')) ~ In" ' 62 |x - x'|. (53) 

5.3.4. N > 3 

Numerical analysis of Eq. (^7|) shows that solutions, compatible with Eq. (pl7|), 
are absent at N > 10. Hence, QLRO is absent for any N > 10. This agrees with 
the previous results for the spherical model.0'i~i For each integer N < 10 the 
RG equation (^) has exactly one solution satisfying the inequality (47). These 
solutions are described in Table 1. In the table, r\ is the critical exponent of the 
two-spin correlation function, the scaling dimension of the temperature (pi). 

Unfortunately, it is not clear if the fixed points, found at N > 3, survive in 3 
dimensions. A zero-temperature fixed point can exist only if the scaling dimension 
of the temperature is negative. Table |l| shows that scaling dimension is positive in 
the one-loop approximation at e = 1 and N > 5. In the 3-dimensional 0(4) model 
the one-loop correction to the scaling dimension —2(N — 2)i?"(0) w 0.7e is close 
to the zero-loop approximation —2 + e. Thus, the next orders of the perturbation 
theory are crucial to understand what happens in 3 dimensions. 

In the 0(2) model the scaling dimension At = — 2 + e is exact. E3'c3 Hence, 
QLRO disappears in 2 dimensions. In systems with larger numbers of magnetization 
components fluctuations become stronger. Thus, one expects the absence of QLRO 
in all the two-dimensional O(N) models. 

At zero temperature Eq. (|27| ) is valid independently of the scaling dimension 
At- It is tempting to assume that at zero temperature, QLRO still exists in the 
RA 0(N > 3) models below the critical dimension, in which At = 0. However, 
the experience of the two-dimensional RF XY model does not support such an 
expectation. Recent numerical simulations show that QLRO is absent even in the 
ground state of that model.L.3 



Table 1. Critical exponents of the RA O(N) model. 



N 2 


3 


4 


5 


6 


7 


8 


9 


r] 7r^e/36 


0.62e 


Lie 


1.7e 


2.7e 


4.6e 


9.0e 


33e 


A T -2 + e 


-2 + 1.3e 


-2 + 1.7e 


-2 + 2.3e 


-2 + 3.2e 


-2 + 4.8e 


-2 + 8.7e 


-2 + 30e 



5.4. Topological defects 

Our RG procedure is based on the decomposition (^o|) which makes sense only if the 
magnetization change is slow at the microscopic scale a. This condition is broken 
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in the core of a topological defect. The nature of defects depends on the particular 
model but we consider different types of defects in the same way.Ej In the XY model 
the defects are dislocation loops. In the Heisenberg model they are point defects. 
Two types of defects are possible in nematic liquid crystals: disclination loops and 
point defects B3 A disclination loop is a line the rotation by 27r around which reverses 
the director: n — > — n. The structure of a point defect is analogous to the structure 
of a hedgehog in the Heisenberg model. 

Topological defects are irrelevant at small e = 4 — D for weak disorder. This 
can be understood from the consideration of the contribution of disclination (dis- 
location) loops and pairs of point defects of size I 3> a to the RG equations at the 
scale I. After averaging over the small-scale fluctuations the size I of topological 
excitations plays the role of the ultra-violet cut-off. The renormalized temperature 
is small: T(l) <C 1. Hence, thermal fluctuations are irrelevant. In the random-field 
problem the disorder-induced term J?(n a n(,) ~ e in the renormalized replica Hamil- 
tonian ( |To| ) is of order (/i 2 ), where the random vector h describes the (renormalized) 
random field E^ 8 = hn, (...) denotes the average over the realizations of disorder. 
In the RA problem the disorder-induced term i?(n a nf,) ~ e is of order (h 4 ), where 
the random vector h describes the (renormalized) random anisotropy E'dis = (hn) 2 . 
Inside a defect the director change is of order 1 at the cut-off scale. Hence, the elas- 
tic excitation energy determined by the renormalized Hamiltonian H (I) ( |l3| ) can be 
compensated by the interaction with disorder only in the positions where h ~ 1. 
The concentration of such positions is exponentially small ~ exp(— 1/e). Thus, 
defects produce corrections of order exp(— 1/e) to the RG equations and do not 
modify the results of the paper qualitatively. The concentration of the topological 
excitations of size I is not more than of order l~ D exp(— 1/e). The above discussion 
is valid, if disorder is weak. In the case of strong disorder, topological defects are 
present at the microscopic scale a and QLRO is absent. Thus, topological defects 
may drive the system into another glass state in which the orientation of the di- 
rector is determined only by the local random potential. The critical strength of 
disorder at which QLRO disappears can be estimated by comparison of the elastic 
and random contributions to the energy. 

The irrelevance of topological defects for weak disorder can also be demonstrated 
with the energy argumentiij modified to take into account the scale dependence of 

fractal structure of large dislocation loops can lead to 



the interaction. A possi 
their strong suppression. 

Recentlythc role of defects in the RF XY model was a subject of intensive 
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discussions. For a review see Ref.t 



6. Comparison with Experiments and Numerical Experiments 

In this section we discuss experimental and numerical results about QLRO in disor- 
dered superconductors, amorphous magnets and nematics in random porous media. 
In the first subsection we consider superconductors. Our discussion is very brief. 
More details can be found in Refs.ErQ The second subsection is devoted to amor- 
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phous magnets. In the last subsection we consider nematics. Since nematics are 
described by the model (§j4|) which differs from the random-anisotropy Heisenberg 
model, we have to derive RG equations for all three Frank constants ifj. It turns out 
that at large scales the difference from the RA Heisenberg model is irrelevant and 
QLRO phase exists at low temperatures and weak disorder. We then make predic- 
tions about light scattering in the QLRO state and discuss several new phases that 
emerge in the presence of external electric or magnetic fields, or mechanical stresses. 
At the end of the subsection we consider existing experimental and numerical data 
concerning QLRO in disordered nematics. 

6.1. Disordered superconductors 

There is a huge amount of literature about vortex states of disordered superconduc- 
tors. For example, in the latest issue of Nature before this review was submitted two 
important experimental observations were published: inverse melting of the voxtex 
lattice in BSCCOtS and discovery of a second vortex liquid phase in YBCO.E3 A 
particular question of QLRO in disordered superconductors has also received much 
attention. 

The most important theoretical predictionEl is the emergence of Bragg peaks in 
scattering experiments with disordered superconductors at weak applied magnetic 
fields. The vortex density p(x) is given by Eq. In scattering experiments the 
square of the modulus \pk\ 2 of the Fourier transform of the density is measured^ 
For a perfect lattice the Fourier transform is non-zero only for k equal to reciprocal 
vectors of the lattice. For a non-ideal lattice we get 

\Pk\ 2 ~ ex P( ik I R i - R/]) exp(ik[iii - ur]). (54) 

ij 

After disorder and thermal averaging one sees that the scattering measures the 
Fourier transform of the correlation function 

C k (x) = ( exp(zk[u(x) - u(0)])) = ( cos(k[u(x) - u(0)])>. (55) 



In terms of the Hamiltonian ( |12| ) for the vortex displacements from the reg- 
ular positions in the Abrikosov lattice, the correlation function ( p9|) is Cj^x) for 
a particular value of k. We have seen that the correlation function (|2^) obeys a 
power dependence on the distance. The generalization of this property for an ar- 
bitrary k is straightforward.tm In terms of the scattering experiment this leads to 
the situation which is intermediate between an ideal lattice and short-range order: 
There areBiagg peaks but they have finite widths. These Bragg peaks were indeed 
observedEjS in agreement with theory. 

Numerical experiments_also support QLRO in the RF XY model of disordered 
superconductors (see Ref£3). 

Finally, we note that the RF XY model is a simplification. A more rea , 
model that takes into account the triangular symmetry of the Abrikosov lattice 
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gives slightly different predictions. An interesting point is nonuniversal QLRO: 
critical exponents are different in different points of the phase diagram. However, 
they vary in a very limited range.HE3 

6.2. Amorphous magnets 

A consequence of QLRO, predicted in the RA Heisenberg model, is divergence of the 
magnetic susceptibility (^). Earlier analysis of Arrott plots, which show the field 
dependence of themagnetization, suggested that the susceptibility does diverge at 
low applied fields.EZI Later it become clear that in strongly disordered amorphous 
magnets the zero-field susceptibility is actually finite (e.g. Ref.Ea). What happens 
at weak disorder is a more difficult experimental question. There is an evidence of 
a finite susceptibility^] at weak random anisotropy. This could be interpreted as an 
experimental argument in favor of the absence of QLRO in the RA Heisenberg model 
in 4 — e dimensions at e = 1. However, in such case the Imry-Ma argument^ and 
the RG analysis in the spirit of section 3 would predict the following scaling for the 
susceptibility: x ~ (J /A) 4 , where J is the exchange strength, A is the anisotropy. 
This scaling was not observed and in contrast to the theoretical expectations it 
turned out that \ ~ * const as A — ► 0. The authors of Ref.S interpreted this as an 
effect of dipole forces. However, the magnetic susceptibility of the pure Heisenberg 
ferromagnet is expected to be infinite even in the presence of dipole forces.O This 
suggests that the experimental system used in Ref.ES cannot be described as a 
realization of the RA Heisenberg model even with dipole forces. Besides, the effect 
of dipole forces on the existence of a QLRO state in the RA Heisenberg model is 
an open question. 

From the numerical side the existence of QLRO is supported by a recent paper El 
However, numerical resultalij support also QLRO in the RF Heisenberg model in 
contradiction with our predictions. We believe that this is a finite-size effect. More 
work is needed for a better numerical understanding of the RA and RF Heisenberg 
models. 

6.3. Disordered nematics 

In the one-constant approximation K± = K2 = K3 the energy F = F$ + F pm 
(||) reduces to the Hamiltonian of the RA Heisenberg model. Since that model 
possesses QLRO the same ordering is expected for the randomly confined nematic. 
However, in all nematics K±, K3 > K2 and this could change the critical exponents 
of the correlation functions in the QLRO state in comparison with the random 
Heisenberg model. Below we demonstrate that this is not the case, i.e. the nematic 
in the porous matrix belongs to the universality class of the RA Heisenberg model£3 
To get a simple idea why it occurs we first consider a two-dimensional nematic film 
with the director n = (n x ,n y ,n z ) = (cos <fi, sin <j), 0) confined in the plane xy of the 
film in the absence of disorder. The Frank energy is Fd = {K\ + -K3XV '<j>) 2 /2 + 
(K3 — ifi){cos 2<fi[(d x (j)) 2 — (d y (j)) 2 ]/2 + sin 2(j>d x <fid y (j)}. The low-temperature phase 
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of this system possesses QLRO, only the term (Ki + K^)(V4i) 2 /2 being relevant at 
large scales since (sin 2(f)) — (cos 2<fi) = 

The systematic consideration is based on the RG equations in 4 — e dimensions. 
Our method follows the line of section 4. All relevant operators of the appropriate 
symmetry are included in the following replica Hamiltonian 



H R = f d 3 r[^ ^(\id a r%d a r% + X 2 d a n a a d^ 

+X 3 nld a n a n a y d^)-J2^^-l (56) 



ab 



where a,b are replica indices, a, /3 = x,y,z label the spatial coordinates, Ai = 
K<z, X 2 = K\ — K 2 , A3 = K3 — K 2 , T is the temperature, the function R(z) describes 
disorder, and the summation over the repeated indices a and j3 is assumed. Due 
to the symmetry n° «-> n° the function R(z) is even. Below we measure the 
temperature in units of K 2, and hence set Ai = 1. To define the energy in 4 — e 
dimensions we add to the Hamiltonian (|56| ) the term Ao ^2 a p d a Ti%d a ng/2, where a 
labels the coordinates in the (1 — e)-dimensional subspace, /3 = x,y, z. The stability 
conditions-^ K\ , K 3 > lead to the inequality 



A2, A3 > 



(57) 



At each step of the RG procedure which is exactly the same as in section 4 
we require that Ai = 1 is unchanged. We get two additional RG equations in 
comparison with the RA and RF Heisenberg models. These equations describe the 
renormalization of the elastic constants A2 and A3. The RG equations in the first 
order in e = 4 — D read 



where the constant 



dT 
d\nL 



d\ 2 
d\nL 
dX 3 



d\nL 



-(D-2)T+(l-X 3 )C <t> T; 
= -A 2 (l + A 3 )Q ) ; 

-(3A3 + A3 — X2)C < p, 



dR{z = l)/dz 



8vrVAo(l + A 3 ) 



1 + A 2 



describes the fluctuations of the small-scale fields (20) 



(58) 



(59) 



> = (^) = ^ln(L/a). 



(60) 



We omit the RG equations for Ao and R(z) since their structure is irrelevant below. 
Eqs. ( p8| ) have the only fixed point compatible with the stability conditions (|57j). In 
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this fixed point T = A2 = A3 = and Eq. (|5^) reduces to the Hamiltonian of the RA 
Heisenberg model which thus describes the large-distance physics of the randomly 
confined nematic. Since that model possesses QLRO in its low-temperature phase 
for weak disorder, QLRO is also possible in confined nematics. For strong disorder 
or high temperature the ordering disappears. Thus, there are three phases: the 
high-temperature isotropic phase and two low-temperature glass phases with and 
without QLRO. In both glass phases the local orientation of the director is fixed by 
the random potential. The disorder driven transition between the glass phases is 
related with topological defects. 

The large-scale correlations of the director lead to strong small-angle light scat- 
tering. We determine its intensity in the limit of the weak optical anisotropy, i.e. 
assuming that in the dielectric tensor e a/ s = e±S a p + e a n a na the anisotropic term 
e a < ti. In this case the scattering cross-section can be found with the Born ap- 
proximation. The scattering cross-section with the change of the wave vector by q 
is given by the expressionLJ 

<r(q) = \u 2 /(A7rc 2 )i a e a0 (q)f b \ 2 , (61) 

where lu is the light frequency, i and f are the unit vectors specifying the initial and 
final polarizations, e a( a(q) is the Fourier transform of the dielectric tensor. Hence, 
(T (q) ~ (Qq/3( < 3)Qq/3(~ ( i))j where Q a/ 3 = n a np — S a p/3 is the order parameter 
and the angular brackets denote the disorder and thermal averages. In contrast to 
the bulk nematic the scattering is caused not by the thermal fluctuations but by 
the frozen configuration of the director. The cross-section c(q) is proportional to 
the Fourier transform of the correlation function G(r) = (Q Q ^(0)Q a( a(r)). In the 
QLRO state this correlator obeys a power dependence on the distance G(r) ~ r~ r > . 
To calculate the exponent rj we decompose Q a p into small-scale and large-scale 
parts with Eq. ( p0| ) and average over the small-scale fluctuations with Eq. ( p(i[ ) : 

(Q a/3 (O)Q Q/3 (r)) = {<(0)^(0)(1 - 

i 

ij i 

+ 5>a( r )4( r X&&> - W3} = Q' Q/3 (0)Q^(r)[l - 6QlnL/a], (62) 

ij 

where Q' a g = n' a n'^ — Sap/3, (...) denotes the average over the fluctuations of cf>, 
and the relation {4>i4>j) ~ Sij which is valid in the RA Heisenberg fixed point is 
used. The constant — —2R"(0) = 0.309e Eq. (|6p| ) is the same as in the fixed 
point of the RA Heisenberg model (see Eq. (fl9|)). The exponent 77 can be found 
with the iterative use of Eq. ([32|) at each RG step until the scale L = r is reached. 
At the scale r the values of the renormalized director field n' are the same at the 
points and r. Hence, Q' a 8(0)Q' a g(r) ~ 1 and ~ (1 — 6C*0 In L/a) K , where 
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K = ln(r/a)/ln(L/a) is the number of the RG steps. Thus, r\ = 6C0. The small- 
angle scattering cross-section is given by the expressiontj 

<r(q)-q- D+v = <T 4+2 - 9e . (63) 

The uniaxial stress modifies the large-distance behavior. The compression along 
the z-axis can be described by adding to the Hamiltonian the term F$ = An 2 , 
where A > 0, since the deformation tends to make the pore surfaces parallel to the 
xy plane and hence favors the planar configuration of the director. The uniaxial 
stretch is described by F$ — An 2 z with a negative A. In both cases A is proportional 
to the deformation. The effect of the electric field is analogous to the effect of the 
stress but the sign of the electric energyEJ F e = — e a (nE) 2 /87r is fixed for a given 
substance. The RG flow is unstable with respect to the perturbation F s and new 
regimes emerge at the scale R = R c at which the renormalized A{R) ~ 1. The 
critical length R c can be found analogously to the correlation length of the RA 
Heisenberg model in the uniform magnetic field (section 5). At small A the result 
is R c ~ I ^4 1 1 / C 2 2 ) = |^4|-o.5-o.i5e rpj^ s t re t c hed system is long-range-ordered 
at the scales R > R c . The nematic order parameter can be calculated analogously 
to the magnetization of theJBA Heisenberg model in the uniform magnetic field 
and is given by the formul£0 Q = (n a n p - 5 a0 /3) ~ R^ 4, ~ \A\ 0A6 \ Long 
range order can also be achieved by applying an arbitrarily weak external magnetic 
field to the confined nematic since the magnetic contribution to the energytS F m = 
— Xa(nH) 2 /2 has the same structure as the energy related with the uniaxial stretch. 
A more interesting situation emerges under the compression. The director averaged 
over a scale R > R c is confined in the xy-plane. The system is thus described 
by the RA XY model. It possesses QLRO but the critical exponents are different 
from the exponents of the Heisenberg model.lia Thus, at the scale R c the cross-over 
from one QLRO state to another occurs. Using the RA XY fixed point found in 
section 5 and repeating the derivation of Eq. (J33T) one finds the Born light-scattering 
cross-section for q < l/R c - f(q) ~ q~ 4 + e ( 1 + 7r I n the RA XY regime the cross- 
section is anisotropic: The small-angle scattering is suppressed, if the incident or 
scattered light is polarized along the compression direction. The phase diagram in 
the presence of a uniaxial deformation is shown in Fig. 10 

The only experimental attempt to test the existence of QLRoS did not pro- 
vide any decisive evidence. While a power decay of correlation functions was not 
observed, only distances less than the experimentally determined Larkin length (sec- 
tion 3) were probedO However, as discussed in section 3, behavior at the scales, 



which are less than the Larkin length, is the same both with and withoi 



QLRO. 



Besides, the experimentB was done with nematic confined in aerosilJ^ This sys- 
tem cannot be described by the model (|)jj) since disorder is only partially quenched 
and elasticity- mediated nonlocal interactions can exist. The same problem concerns 
a possible application of our results to nematic elastomersJlZEa 

One more problem is that no way to determine the disorder strength in experi- 
mental systems was suggested. In Ref.LJ it was claimed that the effective disorder 
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A/T 

Fig. 1. The phase diagram in the one-constant approximation K = K\ = K2 = K3. h is the 
disorder strength, T the temperature, A the anisotropy. Isotropic QLRO exists in region OXQ 
in the A = plane. Surface OQQi indicates the border between the region of short range order 
(SRO) and the region with anisotropic QLRO. At A < the system possesses long range order 
(LRO). 



strength is less in low-density porous media. This claim however ignores the depen- 
dence of the disorder strength on the structure of the porous media. For example, 
it is evident that no ordering is possible in the case of disconnected pores even if 
the density of the porous media is arbitrarily low. 

Results of numerical experiments are controversial: Ref.E3 supports QLRO but 
RcfEJ does not. We believe that one should take the existing numericaLresults 
with care. In particular, there is a question about finite-size effects in Ref.E3 while 
in Ref£j disorder is not weak. Randomness is introduced only in a small fraction 
of sites, but the random anisotropy on each site is assumed to be infinitely strong. 
This is equivalent to the application of a still-not- weak random field in all six neigh- 
boring sites, and the random fields in these six sites are correlated. The latter also 
makes effective disorder stronger. Hence, the effective disorder strength in the cor- 
responding continuous model is not small. Thus, more experimental and numerical 
work is needed. 

7. Conclusions 

In conclusion, we have demonstrated that the random-field and random-anisotropy 
0(N) models possess quasi-long range order at low temperatures and weak disorder. 
In the random-field O(N) model, quasi-long range order is possible at N — 2. In the 
random-anisotropy model, QLRO exists at N < 10. These results can be applied 
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to vortex states of disordered superconductors, amorphous magnets and nematic 
liquid crystals in random porous media. 

There are many other related systems which we did not consider in this paper. 
As the only examples we mention weakly disordered two-d i m en si on a] XY ferro- 
magnets with both exchange and dipole-dipole interactionsElo and charge density 
waves with long-range Coulomb interactions in two dimensions in the presence of 
disorder £3 These systems can be described as super-quasi-long-range-ordered: The 
correlation functions depend on the distance logarithmically slow similar to the 
correlation function of the four-dimensional random- anisotropy Heisenberg model 

©■ 

The results about glass states of weakly disordered systems of continuous sym- 
metry are a part of a wider development. We have good understanding of weakly 
disordered systems, if disorder breaks only the translational symmetry. Random 
fields and random anisotropy also break the symmetry with respect to transfor- 
mations of the order parameter. In this case much has still to be investigated. 
Even the simplest problem, ran Ho rn -fi eld Ising model,E3 is far from being com- 
pletely understood. Recent resultsEj'E3E3 make it hopeful that a satisfactory theory 
of phase transitions in the random-field Ising model will be developed soon. Such 
theory should include new insights in comparison with the existing theory of phase 
transitions in homogeneous systems and may open important new perspectives in 
statistical mechanics in the presence of disorder. 
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Appendix A 

In this appendix we derive an inequality for the correlation functions of the disor- 
dered systems. We consider a system with the Hamiltonian 

H = y"dx D [ff 1 ( ( Hx))-Mx)m(<Kx))], (A.l) 

where (f> is the order parameter, h the random field with short range correlations, H\ 
may depend on some other random fields. We prove the inequality for the Fourier 
components of the field m: 



G C on(q) < const VGdi^q), (A.2) 



where G d i s (q) = (m(q)m(-q)), G con (q) = (m(q)m(-q)) - (m(q))(m(-q)>, the 
angular brackets denote the thermal averaging, the bar denotes the disorder aver- 
aging. 
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Let P(h) be the distribution function of the field h. Then 



s»W = J ( p C 1 >5j^ m q(' 1 )) D W = 

-J p < h i^^) D ^ (A - 3) 

where J D{h} denotes the integration over the realizations of the random field, 
mq(ft) = / £>•!>} exp(-F/T)m(q) / f £>{>} exp(-i?/T). Applying the Cauchy- 
Bunyakovsky inequality to Eq. ( |A.3[ ) one gets Eq. ( |A.2| ) where the const = 

maxq J] D{h}\dP(h)/dh{q)\ 2 /P{h). 

For systems in the critical domain there is a simple way to understand why the 
inequality is valid not only in the Gaussian case but also in a general situation. 
This is just a consequence of the universality. 
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